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Abstract 

In this paper we find the first and second order perturbations of the induced metric 
and the extrinsic curvature of a non-degenerate hypersurface S in a spacetime (A4,g), 
when the metric g is perturbed arbitrarily to second order and the hypersurface itself 
is allowed to change perturbatively (i.e. to move within spacetime) also to second 
order. The results are fully general and hold in arbitrary dimensions and signature. 
An application of these results for the perturbed matching theory between spacetimes 
is presented. 
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1 Introduction 



The aim of this paper is to analyze how the induced metric and extrinsic curvature of 
a hypersurface change when the spacetime metric is perturbed to second order and the 
hypersurface itself is deformed perturbatively to second order. The motivation to carry out 
this calculation is twofold. Firstly, to study the matching conditions between two spacetimes 
when the metric is perturbed arbitrarily to second order. Secondly, to study the dynamics 
of first and second order perturbations of an (n — l)-brane and its backreaction on the 
bulk. These two problems involve in an essential way knowing the first and second order 
perturbations of the induced metric and extrinsic curvature of the hypersurface. 

A natural question that should be addressed to start off is whether going to second order 
in perturbation theory is necessary and interesting. Taking for granted that perturbation 
theory is useful and powerful for a wealth problems, the point is why second order. First of all, 
our present understanding of linear perturbations (methods involved, subtleties, applications, 
etc.) has reached a level of maturity that allows us to go to second order as a natural 
next step. More importantly, there are many situations where linear theory is not accurate 
enough and second order non-linear effects have to be considered. One fundamental example 
is studying inhomogeneities in our Universe (see [T], 0, [H], jlj, and references therein) 
and particularly structure formation (e.g. [B|,[7|). Another example is anisotropics in the 
cosmic microwave background, where the present and planned sensitives of the detectors is 
high enough so that second order effects can already be detected (or ruled out) (see 0,10], 
[TO] for some recent references). 

Besides cosmology, second order perturbations are also fundamental in slowly rotating 
stars. Since the seminal work of Hartle fl] we know that stationary and axially symmetric 
(rotating) perturbations of static stars have a linear component which sets the star into 
rotation without modifying the shape of its boundary. Its effect to second order is to modify 
the shape of the star as well as to modify the rest of the metric components. Dynamical 
(non-stationary) perturbations of stars (both static or slowly rotating) also require going to 
second order (see ^2] and references therein). Second order perturbations have also been 
applied to back holes, specially to study the close limit in the collision of back holes ([T3]. 
[Tlj and therein references). When compared with numerical work, the results of these 
calculations show that the range of validity of the perturbative regime is much broader than 
expected. 

In brane scenario, perturbations are also relevant. One important aspect has to do with 
stability of branes, for which a proper understanding of how the intrinsic geometry of the 
branes behaves under perturbations is required. Small (first order) perturbations of a world- 
sheet evolving in a fixed flat background were studied in [T5] . The generalization to curved 
background (still fixed) was done in jTH], [IZj- Including the backreaction of the bulk (i.e. 
first order perturbations of a brane in a perturbed bulk) was analyzed in where in 
particular the linear perturbation of the induced metric and second fundamental forms of an 
arbitrary (non-null) brane were calculated. The calculations were carried out using exclu- 
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sively spacetime tensors. This approach is very efficient because there is no need to use two 
different sets of indices (spacetime indices and intrinsic indices on the brane). Linear per- 
turbations of branes have also been studied in the context of cosmological perturbations on 
the brane (see ^H] for reviews). Different approaches have been considered: the most widely 
used involves calculations in explicit coordinate systems (see j2Hj and references therein), 
another approach uses a covariant formalism [23 and a third method is the doubly gauge in- 
variant formalism developed by Mukohyama [22] • The latter is geometrically very clear and 
has the advantage that the gauge freedom in the bulk and on the brane are neatly separated. 
As far as I know [22] and [IB] are the only papers where general first order perturbations of 
an arbitrary non-null hypersurface in an arbitrary spacetime are studied. 

Second order calculations have also been considered for some highly symmetric bulks 
and branes [23] , [21] , [213 , [2H] • So far, the only completely general calculation to second order 
can be found in [2Zj where the Lagrangian density of the Dirac-Goto-Nambu brane was 
calculated when both the bulk and of the worldsheet were perturbed to second order. The 
main motivation in [27] was to find a quadratic Lagrangian for the first order perturbation 
equations. The results in the present paper generalize those in [2Z] and extend the calculation 
of the induced volume form (i.e. the Lagrangian) to the full first and second fundamental 
forms. However, only codimension one branes are considered here, unlike [21] where any 
co dimension was allowed. 

Gauge invariance is a subtle point in any geometrical perturbation theory. In the case 
of hypersurfaces the complexity increases due to the presence of a moving boundary, and 
this has led to some confusion in the literature. Going to second order complicates matters 
both conceptually and in the calculations. Thus, it becomes necessary to state clearly the 
theoretical framework defining the perturbations. In this paper I will use a geometrical 
method that, in some sense, combines the approaches of jTH] and [22] • As in [22], the 
hypersurface and its perturbations will be defined by embeddings of an abstract manifold 
into the spacetime (thus splitting the intrinsic objects from the ambient spacetime objects 
in a clear way). This allows us to separate the two gauge freedom sources neatly. For 
the explicit calculation of perturbations, it is much more economic to use the spacetime 
formalism, as in ^H]. The result is finally written down in terms of intrinsic tensors on the 
unperturbed hypersurface, where all objects are naturally defined. The calculations leading 
to second order perturbations are difficult. If a correct approach is not taken, the expressions 
easily become very large and unmanageable. Thus, I will spend some effort explaining how 
the calculations are performed. A basic ingredient is a Lemma which relates perturbations of 
intrinsic objects to a hypersurface and perturbations of suitable objects in the ambient space. 
This Lemma will be used not only to calculate the perturbations within the hypersurface but 
also to analyze the second order gauge freedom. This approach to the gauge transformations 
is complementary to that in [2B] , [2H] where the so-called knight diffeomorphisms were used. 

The calculations presented in this paper have many potential applications. One which I 
consider very relevant is the matching of spacetimes. Constructing spacetimes from the junc- 
tion of two regions across their boundary has been a very useful tool in gravitational theory. 
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The set of conditions ensuring that two spacetime regions can be joined are well understood 
(see e.g. [3U] for a detailed account). The role of the discontinuity of the second funda- 
mental form and its relationship with distributional parts in the energy-momentum tensor 
were clarified by Israel and have become known as Israel matching conditions [31]. Since 
perturbation theory has been useful for many problems, it is not surprising that perturbing 
spacetimes constructed from matching two regions is also of interest. Obvious examples are 
perturbed stars, voids in the universe (as in the Einstein-Straus cheese model [32]), shells 
of matter, impulsive waves, etc. In this paper I will present, as a by product, the first and 
second order perturbed matching conditions for non-null hypersurfaces. 

The paper is intended to be self-contained, so that all the subtleties and difficulties of 
second order calculations on hypersurfaces can be properly understood. I will, of course, give 
credit to previous results whenever necessary. The paper is organized as follows. In section 
2, a brief summary of perturbation theory is given. This fixes our framework and notation. 
In section 3, the Lemma mentioned above will be stated and proven. This result will be used 
throughout the paper. A consequence of this Lemma is that second order perturbations of 
hypersurfaces can be described by two vector fields defined on the unperturbed hypersurface. 
Their explicit form in an arbitrary coordinate system is discussed in Section 4. Section 5 
deals with first order perturbations of the fundamental forms. Here the results of [18 a and 
[22] mentioned above will be recovered. In this section some useful Lemmas to carry out 
the second order calculations are presented. The result of the second order perturbations 
is stated in Section 6, leaving the details of the proof to Appendices A and B. Section 7 
discusses first and second order gauge transformations. For hypersurfaces there are two 
types of gauge freedom, namely the one coming from the ambient spacetime and another 
one intrinsic to the hypersurface. Both are discussed in this section. Section 8 applies the 
previous results to the perturbed matching conditions between two spacetimes. A theorem 
giving the necessary and sufficient conditions for second order perturbations to match across 
a matching hypersurface in the background is presented. This theorem can be potentially 
applied to many situations. One case that has already been investigated involves first and 
second order stationary and axially symmetric perturbations of spherical stars |33]- The two 
Appendices contain the main steps in the calculations of the second order perturbations. 
The reader who is not interested in detailed calculations may skip the Appendixes and 
concentrate on the main text. 

2 Summary of perturbation theory 

Perturbation theory deals with one parameter families of spacetimes 1 (Ai t ,g t ) and their 
first and higher order variations around one element of the family, say (A^ofl'o)- I n order 

1 A C m spacetime is a Hausdorff, connected C m+1 manifold of dimension n endowed with a Lorentzian 
metric of class C m . If we are considering a manifold- with-boundary then the boundary dM. is also assumed 
to be C m+1 . Our signature and sign conventions for the Riemann and Ricci tensors follow 34 . 
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to take derivatives with respect to e we need some fixed set of points, i.e. we need all A4 e 
to be diffeomorphic to A4 . Through these diffeomorphisms we can define a one-parameter 
family of metrics g e on A4 Q associated to g e . We denote Aio simply as M.. Since we want to 
take second variations of extrinsic curvatures we take A4 to be C 4 , as a manifold. We thus 
consider a differentiable family of Lorentzian metrics on Ai, i.e. a C 2 map 

T : I — ► G(M) 

e — > 9e 

where &(Ai) denotes the set of C 3 symmetric, non-degenerate, two- index tensor fields in M. 

can be endowed with a natural differential structure, T is C 2 with respect to this 
structure). 

We denote by V e the Levi-Civita covariant derivative of the metric g e . Given two arbitrary 
metrics g and ~g, the corresponding covariant derivatives V and V are well-known to be 
related by 

V,T; = V M T« + C^T: - C^, (1) 

where C^ 7 = (l/2)g a/i (V pg^ + V 7 ^ M/3 — V^g^). Similar expressions hold for tensors with 
any number of indices. Substituting g by g t= o and g by g e , we can take e-derivatives in (JTj) 
at e = to get 

O/f rpa rpV fry\ 

— ° fry, v — ° u^-tpi l Z J 

6=0 

where we have defined 

, = \ (V^ + V 7 i^ - V«K' M ) . (3) 

Here and in the following we set go — > g and V° — > V. The tensor K' is the first order 
perturbation of g along the family g e . We will simply call it first order perturbed metric. 
Similar expressions as (0) hold for tensors with any number of indices. The second derivatives 
of (JTJ) gives 

= + T? - (S"l + 2K' a s S%) Tj£, 

e=0 

where K" a p = d |^ | £= o is the second order perturbed metric and 

S'% = \ (V^"° 7 + V 7 K% - V Q n 7 ) . (4) 



a/3 



<jge a/: 



de 2 
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3 A useful Lemma 



We want to calculate how the first and second fundamental forms of a hypersurface change 
when the ambient metric and the hypersurface are perturbed. Thus, we consider a one- 
parameter family of hypersurfaces E e of A4. As before, in order to define variations we need 
a fixed set of points, so we assume all E e to be diffeomorphic to each other. We allow the 
hypersurfaces to change both as a subset of points of A4 and also in the way we coordinate 
them for different e. Since the fundamental forms are pull-backs of tensors on A4, their 
dependence on e arise because of three facts: (i) because the ambient metric depends on e, 
(ii) because the hypersurface E e considered as a subset of A4 changes with e and (iii) because 
the way in which we coordinate E e is allowed to depend on e (even if the hypersurface as a 
set of points remains unchanged). Points (ii) and (iii) can be treated together by viewing the 
hypersurfaces E e as embedded hypersurfaces, i.e. as the images of a family of embeddings 
$ e : E — > A4, where E is a copy of any of the E e , say E . It is useful to view E as an abstract 
manifold detached from the spacetime so that one knows clearly where the different objects 
are defined. Thus, we shall distinguish between E as an (n — l)-dimensional manifold and 
E = $o(S), which is the hypersurface in A4. 

The fact that the fundamental forms on E depend on e from several sources and that we 
want to do the calculation up to second order makes it very important to use a method as 
covariant as possible and use coordinates only when absolutely necessary. Moreover, it is very 
convenient to work on the ambient manifold as much as possible and perform the pull-back 
only at the very end (in agreement with [TH]). The alternative of calculating the derivative 
directly on E is, of course, possible but much more difficult. In this section I present a Lemma 
which shows how derivatives of geometric tensors on an embedded submanifold (or arbitrary 
codimension, including codimension 0) can be calculated from derivatives performed purely 
on the ambient manifold. This result will be crucial for the calculations in the following 
section. Throughout this section all differentiable objects are C 3 unless otherwise specified. 

Thus, let Af and A4 be two differentiable manifolds and let Xe '■ Af — > A4 be a family 
of differentiable maps. Let us consider a C 2 family of covariant tensor fields T e on A4. We 
can pull-back this family to Af and define a one-parameter family of tensors T t on Af. We 
are interested in determining the first and second derivatives of % with respect to e. Using 
directly the definition of derivative, 

= nm xU h (Te + h)~Xt (Te) = Um xU h (Te +h )-xt {T t+h ) ^ / T e+h - T e \ 
de h™o h h™o h Xe \h™o h J 

where we have added and subtracted x*(T e+/l ) in the numerator and we have used the 
linearity of x*- The second term is the pull-back of the derivative of T e . The first term cannot 
be written directly in a simple form because there is no a priori relationship between x*+h anc ^ 
x\ (contrarily to what happens for instance in a one-parameter group of diffeomorphisms) . 
Assume now that Xe are embeddings. Then, there exists a set of diffeomorphisms : A4 — > 
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A4 of the ambient space such that the diagram 

AT 

Xh+e 




(6) 



is commutative. We moreover fix \1/q = 1^ (the identity on A4) for all e. Geometrically, 
^ e h transforms the point Xe(p) into the point Xe+h(p) for all p G Af. Notice that ty e h is fixed 
only on Xe(Af) and therefore it is non- unique in general. However, when M has the same 
dimension as A4 and Xe are diffeomorphisms, then ^> e h is unique and given by Xe+h ° X7 1 - 

The following Lemma gives explicit expressions for the first and second e-derivatives of 
T t in terms of objects defined in the ambient space A4. 

Lemma 1 Let T e be a C 2 one-parameter family of covariant tensor fields on A4, Xe '■ Af — » 
AA a C 2 family of embeddings and define T t = x*(T e ). Then 



d 2 T e 
de 2 



Xe \ <£wT e + £^ c £^T t + 2£ 2 t ( 



fdT e 



\ de 



+ 



d 2 T e 
de 2 



(7) 
(8) 



where 



dm 



dh 



h=0 



dZ e 
de 



(9) 



and $? e h is any set of diffeomorphisms of AA. which makes the diagram (Q|) commutative. 

Proof: The commutativity of the diagram (jUJ) implies that the first term in (jSJ) can be 
written as 



lim 



Xt +h (Te- 



+h 



h 



X'e 



lim 



*t (Te 



+h) 



h 



Xt (£zje) 



where the definition of Lie derivative has been used. This proves (|7j). For the second 
derivative we apply this expression twice 



d 2 % 
de 2 



d 



d 



de J V Ze de r 

r fdTA d(£ 2 T e ) d 2 T € 

Xl , £,£,T e + £,[— + + — 



(10) 



Only the third term needs to be elaborated. Using again the definition of derivative and 
adding and subtracting a suitable term we get 

de h->o h h-*o h 

which, using the linearity of the Lie derivative, becomes 

d (£zj,) _ fdT e 
de ~sr e \ de 

Inserting this into (fTUj) and using the definition W e = ^ the lemma follows. □ 
Remark. Notice that when M is a hypersurface E, all the information regarding the 
first and second variation of the hypersurface E e around E is encoded in the two vector 
fields Z\ = Z e=0 and W = W e=0 . These vectors are, by construction, defined everywhere 
on A4 (because they are defined in terms of the diffeomorphisms ^ e h ). However, only their 
values on So should matter. Geometrically, they define how the hypersurface is deformed to 
first and second order (as we shall see, the second order variation is best defined by using a 
suitable combination of W and Zi). These vectors have tangential and normal components. 
The normal part determines how the hypersurface moves in spacetime as a set of points, 
while the tangential part encodes the information on how the different S e are coordinated. 
The fact that these vectors have been extended off S , and that this extension is essentially 
arbitrary (due to the large freedom in defining ty e h ) provides a powerful check for the validity 
of the final results, namely that they must be independent of the extension of these vectors 
outside S . 



4 First and second order perturbations vectors of E 

Let us now concentrate in the case where M is a hypersurface S. We replace Xt — > &e i n ah 
the expressions above. Our aim in this section is to find explicit expressions for the vectors 



Zi 



an 

dh 



=h=0 



and W = 4g* 

de 



, so that they can be determined in explicit examples. For 



the second variation we use Z 2 = W + ^7 g Z-y instead of W. The reason will become clear 

later on. We call Z\ and Z 2 respectively as first and second order perturbation vectors o/E. 
Let us choose local coordinate systems on E and on M. so that $ e are written 



$ e : E — > M 
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Proposition 1 The first and second order perturbation vectors Zi a (y) and Z^iy) of the 
hypersurface £ read 



Zi a (y l 



d<$> a {y\e) 



de 

d 2 $ a (y\e 



6=0 



de 2 



'ID 



;i2) 



6=0 



where Xo(y l ) is the local form of the unperturbed embedding $ - 

Proof: : M — > M satisfies $ e+h = ^f e h o $ e . Let ^ a (x a } h, e) be the local coordinate 
form of ^ e h . By construction, $ a and satisfy 

^{y\e + h) = ^ a {^{y\ e),h,e), 

which has as immediate consequence 



(13) 



d^ c 



dh 



d^ a 



dx 13 



(3>(yV),V) 



de 



+ 



d^ c 



(l/V) 



de 



($(y\<i),h,e) 



Evaluating at e = h — and using the fact that is the identity on Ai for all e we obtain, 



d^ c 



dh 



d$ a {y\e) 



de 



(x a =xfi(y i ),h=0,6=0) 

According to its definition, Z e has components Zf(x 



e=0 



dh 



h=0 



(14) 



. Expression (JTTj) 



follows directly from ()14p. For the second order perturbation vector Z 2 , let us first find the 
coordinate form of d e Z e + V % Z e . Directly from its definition one finds 



d e Z?(x) + V* Zf{x) 



d 2 ^ a 



e=0 dhde 
d^ 



+ 



0,0,0) 
9^7 



(x,0,0) 



dx^dh 



(as,0,0) 



(as,0,0) 



(15) 



0,0,0) 



which contains no second derivatives with respect to h. Performing the second e, h derivative 
of ((EH) and evaluating at e = h = gives 



d 2 $ a ( y \ e + h) 



dedh 



d 2 ^ a 



e=h=0 



dx^dh 



d&(y\, 



{x=x {y i ),0,0) 



de 



d 2 ^ a 



+ 



£=0 



dhde 



(16) 



(x=x (y'),0,0) 



Taking into account d e dh$ a (y\ e + h)\, e=0 h=0 ^ = d 2 § a (y\ e + h)\, e=0 h=0 \, and noticing that 
d 2 $ a (y\ e + h)\ e=h=0 = d e d e $ a (y\ e)| e=0 , the vector field d € Z« + V gZf evaluated on £ and 
at e = becomes 



d e Z« + V?Z, 



d 2 $ a (y\e) 



(x=x (y i ),e=0) 



de 2 



6=0 



+ r^(x Q ( y i ))z?( y i )zi( y % (17) 



and (H2J) follows directly from its definition Z 2 a (y i ) = d t Z* + ^ 2 € Zf\ □ 

Remark: Both the first and second order perturbation vectors Z\ and Z 2 depend only 
on the family of embeddings $ e (i.e. on the hypersurfaces S e ) and not on the specific 
choice of ^ e h . Moreover, it is clear from (|T2j) that Z 2 a (?/ 1 ) corresponds to the covariant 
acceleration, evaluated on S , of the curve defined by the motion of a fixed point of £ when 
the hypersurface moves, i.e. by the curve Q a (y\e) with y l fixed. Notice also that, had we 
chosen W = d e Z e \ t=0 as our second order perturbation vector, we would have obtained a 
vector field which depends on ty e h , i.e. it would not be defined solely in terms of the one 
parameter family of embeddings. This is why Z 2 is preferable to W. 

5 First order perturbation of the hypersurface 

We can now start the calculations of the first and second order variations of the fundamental 
forms of S . We shall assume that this hypersurface contains no null points, i.e. that its 
induced first fundamental form h = $o(fiO defines a metric. For small enough e the same 
will be true for E e at least on compact subsets. Since we are only interested in derivatives 
at e = we can assume without loss of generality that all S e are non-degenerate. Let us 
denote by h e = &*{g) the one-parameter family of induced metrics. Notice that all of them 
are defined on the same manifold E. Let also 2 n e be the unit normal to S e with respect 
to g e . Its orientation is taken arbitrarily on So an extended to all e by continuity. Let us 
extend n e to an open neighbourhood U of S e . By working locally near one point we can, 
without loss of generality, choose U to be independent of e. We keep n e unit everywhere on 
U, i.e. geapn^n^u = a, where a = +1 for timelike hypersurfaces and a = —1 for spacelike 
ones. Hence, we have at hand a one-parameter family of one-forms n e defined everywhere 
on U and both covariant derivatives at constant e and e-derivatives at fixed spacetime point 
x G U can be performed. £ inherits a one-parameter family of second fundamental forms 
n e = <3>*(V e n e ). We drop the subindex for any background object, thus E is endowed 
with a metric h, covariant derivative D, second fundamental form k and has unit normal 
n. We also write the background embedding simply as All spacetime indices are lowered 
and raised with the background metric g a p and it inverse. Similarly all hypersurface indices 
are lowered and raised with hij and its inverse. Tensors on £ will carry Latin indices. 

The first and second order perturbations of the induced metric and second fundamental 
forms are obviously 

First order perturbations: h' 

Second order perturbations: h h 



2 In this paper boldface letters are used to denote one- forms. 



dh € 
Ik 
d 2 h 



0K € 

Ik 



de 2 



d 2 n f 



K 



de 2 
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These derivatives are taken at fixed point p in the abstract manifold E. In this section we 
shall obtain the explicit expressions for the first order perturbed quantities. The second 
order quantities are considered in the next section. 

Let us start by introducing some notation. It is well known that covariant tensors on the 
background hypersurface E are in one to one correspondence with spacetime tensors defined 
on E and which are totally tangent to E , i.e. tensors C^ 1 ... Atm satisfying C fll ...^ m n pa = (a = 
1 ■ • ■ m) on E . The one- form n is obviously hypersurface orthogonal on E . We can assume 
without loss of generality that its extension off E is chosen so that this property is kept. 
We could also choose n so that it defines a geodesic affinely parametrized congruence. These 
two conditions would fix n uniquely and would simplify the calculations below. However 
we prefer to leave the acceleration a = Van completely free. The increase in complexity is 
compensated by the fact that the result has to be independent of a. This provides a non- 
trivial check for the validity of the result. Since, as we shall see, the calculation is quite 
involved, it is convenient to keep non-trivial checks at hand. 

Being n hypersurface orthogonal, its covariant derivative reads 

V a np = an a a[3 + K af3 , (18) 

where K a p is symmetric and completely orthogonal to n. This tensor is obviously the 
counterpart on Ai of the second fundamental form Kij. From ()18|) it follows £ag a i3 = 
an a af3 + an^a a + 2n a p. 

Covariant derivatives of a tensor C ix ... im within the hypersurface can be calculated by 
considering its counterpart C Ml ... Mm on spacetime. Indeed, if we extend this tensor to a 
neighbourhood of So in such a way that it remains orthogonal to n (and otherwise arbitrar- 
ily), the three dimensional covariant derivative can be calculated as 

D a C^...^ m = h^h^l ■ ■ ■ h^VvCp-,...^, 

where = 5^ — on a np is the projector to the hypersurface. More concretely, this means 
that DjCi v ..i m h as DaC^...^ as its spacetime counterpart. A simple integration by parts 
shows that covariant derivatives and three-dimensional derivatives are related by 

n 

(19) 

i=l 

For later use let us notice some useful expressions. The first one is obvious: for any 
covariant tensor A ai ... am and any function F 

n 

(£FnA) ai ... am = F (£nA) ai ... am + A ai ...,... am n p V ai F (20) 
Less immediate, but still easy, are the following three Lemmas. The first one is well-known 
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Lemma 2 Let £ be an embedded submanifold of Ad with embedding $ : £ — > M. and 
let V be a vector field on a neighbourhood U of tangent to this hypersurface (i.e. 

= ^*(^e) for some vector Vs on E j. Then, for any covariant tensor A onU 

<T (£?A) = £ % ($M) . 

Remark. For simplicity we will use the same symbol to denote Ve and V. The precise 
meaning will be clear from the context. 

A consequence of this Lemma is that, for any vector field X = Rn + V, with V orthogonal 

to n, 

$* (£~ g ) = £~ h + 2Rk, (21) 

Lemma 3 Let B a p be any symmetric tensor and X any vector field. Defining S(B)^ = 
\{V a B^ + VpB^ - WBap) and H a = B^X^. The following identity holds 

(£ j[ B) ap + 2X ll S(B)^ = (£ s g) all . 

The next Lemma tells us how to perform second Lie derivatives. The Riemann tensor of 
(A4, g) is denoted by Ra^s- 

Lemma 4 Let X be an arbitrary vector field and B a p any symmetric tensor. Then 

( £ x £ x B ) Q = ( £\j?x B ) + X^X V V 'pV v B a p - B av X lx X~ f R v . p - B^X^X^ R v 

+2 {X^V^BJ) VpX" + 2 (X^^) V a X» + 2B, V (V Q A^) (V p X v ) , 

Proof: Expand the first term and use the Ricci identity applied to X. □ 
A particular case of this Lemma is obtained for B a p = g a p: 

(£ x £ x g) aP = {£ WgX gU ~ ^x v R^ v + 2 (v a x^) (v^) . (22) 

Combining this with the general expression (|2T)|) it follows easily that, for any pair of functions 
F x and F 2 , 

£ F-i_n£ F 2 n g a f3 = £F 1 F 2 a9al3 + 2F\F 2 (— U^n" R a ^l3u + Kan K /3J 

+a (V a Ft V/jF 2 + V^VpFx) + 2/€ a/3 F 1 n(F 2 ) + n a Gp + n p G a , (23) 

where G a = aa a F 1 n(F 2 ) + £F 1 n{V a F 2 ) + 2aF 1 F 2 a fl K^ a + F^n^a^a^. With F 1 —> Q 1 and 
F 2 — > 1 we get 

£ Q\n,£ n9a/3 £ Qia9a/3 

+ 2Qi (-nVRa/ipv + Ka»Kf? + naUpa^ + an a a^Kp^ + anpa^K^) (24) 
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It is convenient to decompose the first perturbed metric into tangential and normal compo- 
nents with respect to n. Explicitly 

K' a p = Y'n a n/3 + an a r' ' p + anpr' a + K' a/3 , (25) 

where obviously r' a and K' l a p are orthogonal to n a . We can now find the first order 
perturbations of the first and second fundamental forms of E . 

Proposition 2 (Battye & Carter, 1995) Let (M,g) be a C 3 spacetime of any dimension 
and S an arbitrary non-degenerate hypersurface defined by an embedding $ : S — > M.. Let 
h be the induced metric, k the extrinsic curvature and n the unit normal to the hypersurface. 

If the metric g is perturbed to first order with K' and the hypersurface is perturbed to 
first order with a vector field Z\ = Q x n + T X) where T\ is tangent to E , then the induced 
metric and extrinsic curvature are perturbed to first order as 

h' i:j = £fh ij + 2Q l K ij + K' af3 e^e l3 j , (26) 
n'ij = £f i n ij - aDiDjQx + Qi (-n^n" Ra^e^e 1 ^ + k u kJj + ^Y'k^ ~ n nS'aj3 e i e j i 

where Y' = K' al3 n a n 13 , S f is given in J3J) and ef = are tangent vectors to E . 

Proof: From Lemma ^ 

h' = d e h e \ €=0 = (£ z ~ i9 ) + $* (K') (27) 

and (|27)|) follows directly from LemmaEland ()21jl . For k', we notice that 2n e = 2$*(V e n e ) = 
§*(£ri e g e ). Applying Lemma [T] 

2d e K e \ e=0 = <T (£ fi £ H g + £^g + £ n K') , (28) 

where n*i = d e n e \ e=Q . Let us identify this vector: for its normal component we use that 
n e is unit for all e, i.e. (n t ,n e ) g<£ = a. The derivative at e = gives n\ a n a = — ~Y r . 
For its tangential part, it is convenient to use the normal one-form, m x = d 6 n e \ t=0 . From 
$*(n e ) = 0, Lemma[T]gives $>*{£g n + m x ) = 0. Using £ft(n) a = a a and (J2T)|) with A = n, 
F = Qi gives 

V{m x ) i = -{Q x a i + aD i Q x ). (29) 

The identity d e g^ = —g^g^ v d e g e(iu implies that n\ and mi are related by n\ a = —K' a ^np + 
mi a . Hence mi a n a = |Y' which, together with (|2l?|) . gives 

mi a = ^Y'n a - (Q ia a + aD a Q l ) , (30) 
ri! a = ~Y'n a - (r' a + Q x a a + aD a Q 1 ) , (31) 
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where the decomposition has been used. Inserting (pTT|) and (jUj) into (|25j) yields 

29 e /« e y| e= Q = 2£rp^Kij + £Qiahij + 2Ql n U R ap p u e^ C^ + KnKj\ — (jY Kij 

-£? +Ql a + ,^ Ql hij + (i?«Jf) , (32) 

where (gradQi) 1 = D l Q x . For the last term we apply Lemma El with B a p = K' a p and 
X = n, i.e. 

(£ H K') af3 = + £p +aYlfi g a p (33) 

and the Lemma follows. □ 
Remark. This result was also derived by Mukohyama [22] hi his study of first order 
perturbed matching conditions between spacetimes. The proof presented here is based on 
Lemma n which makes the calculations very efficient. This will allow us to push the calcu- 
lation to second order. 

6 Second order perturbation of the hypersurface 

Lemma 0] gives an expression for second order Lie derivatives. The next Lemma gives an 
alternative expression adapted to decompositions into tangential and normal components to 
the hypersurface. 

Lemma 5 Let S be an arbitrary non-generate hypersurface in a spacetime (A4,g). Let n 
be a hypersurface orthogonal unit vector in a neighbourhood o/Eq which is orthogonal to E . 
Let B be any tensor and Z x = Q x n + T x , with T x orthogonal to n. Then 

£ - Z\^ £ 'f 1 £ Ti^ ^£ T 1 £ QinB £ Cin+Di^ £ Q\a,B 4~ £ Q 1 ft£ Q\nB , (34) 
' " 1 z i 

where 

Ci = Qxn{Q x ) + 2T 1 (Q 1 ) - aT x a T x P K aP , = 2Q 1 T 1 a n p + T x a D a T x p . (35) 

Proof: Applying (|19p to Ti„ in the decomposition Z x = Q x n + T X , we get 

V Q Zi M = an a (Qxdp + n p V p T lfl ) + n M (V a Q x - aT x p K pa ) + Q X K afi + D a T lfl . (36) 
Furthermore, linearity and the general property £^£f — £f£^ = £[x?] hnply 

£ ' z 1 £ z x B £ rp^£ rj^B -\- 2 £ j^£ q^hB ~\- £iQ^rpiB -\- £ q 1 ^£ q^B . (37) 

We want to introduce a term ^ g Z x in the right hand side. To that aim we rewrite [Q x n, T x ] 
in the following way 

[Qxn^f = [Z X ,T X }» = VgTf - V f Z x » = V^Zf - V fi (Qm») - v f z x » = 

V^Z/ - Q\a p - 2Q 1 T 1 Q < - T x a D a T^ + ra" (-Q x n{Q x ) - 2T X (Q X ) + ctT^T/k^ , (38) 
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where the first three equalities are immediate and the last one follows directly from (JTSj) and 
(jnHJ). Combining (|S7|) and (jHHj) the Lemma follows. □ 
As before, let us decompose the second order perturbed metric into tangential and normal 
components 

K" aP = Y"n a n p + an a r% + on p r" a + K" l a p. (39) 
We can now write down our main result of this section. The proof is given in Appendix A. 

Proposition 3 With the same assumptions and notation as in Proposition^ if the metric 
is perturbed to second order with K" and the hypersurface is perturbed to second order with 
Z 2 = Qin + T 2 ( with T 2 orthogonal to n) then the induced metric and extrinsic curvature are 
perturbed to second order as 

h ij £ rp^hij ~\~ 2Q 2 Kij ~\~ a{3^i ~t~ 2£j^h ij £rp^£q=^hij -\- 
Jr£ 2Q 1 r'-2QMf 1 )-D T - i f 1 h ij + 2 (o'TiT^Kis - 2T 1 (Q 1 ) + 2aQ{Y'^j Kij + 

+2QI (-nWR^e^ + k u k 1 ^ + 2aD i Q 1 D j Q 1 - AQ^S'^efe^ (40) 

n'ij = £f 2 Kij - aDiDjQ 2 - Q 2 n^n L 'R a ^ v e^e P j + Q 2 KuK l k - n^S^e^e? 
+2£ fi n' i3 + K i:j ( |y" - -Y' 2 - a (n + trAQi) (r l + crD l Q 1 ) + 2aQ 1 n li n p n 6 S 
+ (aY'n^ + 2r M + 2aD fX Q 1 ) - 2Q 1 n li n" '(V ' - 2nX^'LefD J Q 1 

+ ^ CT gr a d(T 1 (Q 1 ))-igrad(T 1 ; T 1 m K im )-|y'g r ad(Q 1 ) + 2 CT Q lK ( gr adQ 1 ) ^ij 

+ (2T 1 (g 1 ) - uT x l T x m K lm - aQxY'^j (nWR^efe* - k u k 1 ^ + I {D i Q 1 D j Y'+ 

-\-DjQ\DjY ) — £f^£f^Kij — £ 2Q 1K (T 1 )+Df T\ K ij ~ 2o r Q 1 <£ grad (Q 1 )Kj J ' 

-Ql (nWn s (\7 5 R a ^u)e?eP + 2nWR Smv efe'*K l i + 2n»n h 'i^efe? 4) , (41) 
where S" is given in and, /or any tangent vector V , k(V) is the vector with components 

Remark: Propositions El and El are still true for metrics g of arbitrary signature. 

The expressions in this Proposition are rather involved, and the calculations leading to 
them are not easy. Thus, it is useful to have ways of checking whether the expressions are 
indeed correct. We have already mentioned two such tests, namely that the result must be 
completely independent of the extension of the normal n and of the perturbations vectors Z\ 
and Z 2 outside Sq. This is obvious from the expressions above because no term containing 
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the acceleration a is present. Moreover, there are no normal derivatives of objects defined 
intrinsically on E . Another important test is that the expressions must transform correctly 
under gauge transformations. In the next section we study gauge transformations when 
there is a hypersurface present. 

7 First and second order gauge transformations 

The definition of the perturbation of the metric and of the hypersurface is based on an a 
priori identification of Ai e with A4o — A4 and of S e with the abstract manifold E. Obviously 
this identification is not unique. The consequence of this non-uniqueness if the gauge freedom 
inherent to any geometric perturbation theory. In our case we have two different sources of 
gauge freedom, namely the identification of the ambient space and the identification of the 
hyper surf aces. Let us describe them in detail. 

The freedom in the identification of A4 e is described by the possibility of performing an 
arbitrary diffeomorphism of M. t (which, of course, will in general depend on e) before doing 
the identification with M. = Ai . An equivalent way of stating this fact is that for any 
one-parameter family of metrics g e on M. we can define an infinite number of equivalent 
families by performing an e-diffeomorphism of g e . The equivalence is shown as follows. 
Each Ai e has a metric that we denote by g e . These metrics cannot be compared with each 
other because they are defined on different spaces. Thus, we need a one-parameter family 
of diffeomorphisms A e : M. — > Ai e in order to be able to relate them, and this defines a 
one-parameter family of metrics on M. by g t = A*(g t ). If we now perform a diffeomorphism 
fi e on Ai before identifying with A4 e , it is clear that A e becomes A g) = A e o n e (the 
superscript (g) stands for "gauge transformed"), and the new parameter family of metrics 
is = A^*(g e ) = Q*(g e ), as claimed. All these families are geometrically equivalent. 
So, although the first and second order perturbation metrics obtained from gi 9 ^ and g e are 
indeed different, they are intrinsically the same (i.e. gauge equivalent). In this section we 
use Lemma [T] to find the gauge transformation law for K' a/3 and K" a/3 (c.f. j2H] and [2H] 
where the so-called knight diffeomorphism are used to describe gauge transformations of any 
order) . 

Let us then apply Lemma ^ with M = M. and Xe = In this case the diffeomorphisms 
which make the diagram (0) commutative are uniquely given by ^ e h = Qh+e ^ 1 - Choose 

T e = g e so that % = gi 9 \ We write, as usual, K' = ^| e=0 , K" = ^^| e=0 and we define their 

gauge transformed tensors as 




dgt 
de 



K" = 

9 




de 2 



e=0 



e=0 
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Let us also define 



Si 



dtt t 
Ik 



h+e 

~dh 



n 



v 



h=0 



dS e 
~ck 



S 2 = V + V*S 1 



(42) 



e=0 



and call S\ and S2 the first and second order gauge vectors. 

Proposition 4 (Bruni et. al., 1997) Under a gauge transformation defined by the vec- 
tors Si and S 2 , the first and second order perturbed metric transform as 



K" 

9 a/3 



K" a /3 + £ s 2 9a/3 + 2£ ^ ^ a/3 ~ 2Si^ Si' ' R a ^f3v + 2 V a 5'i A1 V^i 



Proof: The first expression is a direct consequence of Lemma ^ For the second order 
perturbation, the same Lemma gives 

K" g = £ v g + £ gl £ gl 9 + 2£ g K' + K" . 

The Proposition follows directly from Lemma 0] with X — > Si. □ 
Remark: In |2Hj the vector V was used instead of S2 to define second order gauge 
transformations. Both vectors are equally suited in this case (unlike for hypersurfaces, 
where Z2 is intrinsic to the perturbation E while W is not, see Remark after Proposition HJ. 
Second order gauge transformation were analyzed also in [2^ and S2 was used there. 

When performing a gauge transformation not only the perturbed metrics change but 
also the vectors Zi and Z2 of the hypersurface are modified. Geometrically, this is clear 
because changing the way how the different manifolds Ai e are identified to each other affects 
how the abstract manifold E is embedded into M. at different e, and this obviously changes 
Z\ and Z 2 . It is clear, for instance, that in suitably chosen gauges one can always make 
these two vectors identically zero. Let us therefore determine the behaviour of Z\ and Z 2 

under general gauge transformations. We denote Zi 9 and Z 2 ^ the corresponding gauge 
transformed vectors. 

Proposition 5 Under a gauge transformation defined by Si and S2, the first and second 
order perturbation vectors of E transform as 



(9) 



-(g) 

1 — ^1 — Oli 
Z 2 - S 2 - 2V fi Si + 2V s/i. 



(43) 



Proof: The hypersurface E e is intrinsically (i.e. gauge independently) embedded into M. e . 
However, we are viewing them as hypersurfaces of M. via the embeddings <3> e : E — » M.. 
Recalling that the identification of M. and Ai e is given by A e it follows that A e o $ e is gauge 
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independent, i.e. A^ o $^ = A e o $ e . Under a gauge transformation A e transforms to 
A^ = A e o f2 € , therefore the transformation law for the embeddings is 

= Q- 1 o 



Hence, the gauge transformed diffeomorphisms which make the diagram commutative (jOJ) 
(i.e. such that &l+ h = ^jf e o <E>^) can be chosen to be e = ° ^ ° ^e> which is more 
conveniently written as ^ e h oO e = f2 e+ h o ^jf . Taking the first derivative with respect to h 
at h = we find 



Z e (a(aO) = ^e(x) + 



dx a 



z (a)< 



x 



(44) 



where RJx) 



dO £ Q) 



;h . This vector is related to S t in (J42j) by 



S e (x) 



dh 



R e {n-\x)). 



(45) 



h=0 



(<?) 



Putting e = in (j44j) and recalling that Q Q = I_m the first order transformation Z\ = Z\—S\ 
follows. For the second order, we take an e-derivative of f)44j) at e = 0. Applying directly the 
definitions we find 



W a {x) + S l p {x)dpZ l a {x) 



dR?(x) 



de 



Now, the derivative of (|4o|) at e = gives 



+ Z 1 (9 ^{x)d /3 S 1 a {x) + W {9)a {x). 



(46) 



£=0 



dRf(x) 



d 2 Q a 



e=0 



de 2 



V a (x) + S 1 /3 (x)d S 1 a (x). 



(47) 



e=0 



Inserting this into (|47)|) and recalling the definition of S 2 and Z2 the transformation (|4Hjl 
follows. □ 
Remark: From (|47|) a useful expression for S2 directly in terms of derivatives of Q e 
follows (c.f. Proposition 



S 2 a (x) 



d 2 Q?(x) 



de 2 



e=0 



+ T a ^(x)S 1 P(x)SS(x) 



Under spacetime gauge transformations, the tensors h', h", k', and k" must be gauge 
invariant because they are defined intrinsically on E. This provides a very strong potential 
check for the validity of the expressions given in Proposition El The calculations required 
to perform the check are however very involved and have not been done analytically. Nev- 
ertheless, with the aid of a computer algebra program written in Reduce, I have checked 
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gauge invariance for an important number of non-trivial examples, with positive results in 
all cases. 

As already mentioned perturbed hypersurfaces in a perturbed spacetime have two in- 
dependent gauge transformations. The first one has already been analyzed. The second 
gauge freedom comes from the fact that the hypersurfaces E e embedded into Ai € must be 
identified with an abstract copy of E. This identification entails the freedom of perform- 
ing an e-dependent diffeomorphism Xe of £ before embedding this manifold into M.. Thus, 
the gauge freedom is given by the transformation = $ e o Xe (recall that $ e : £ — > M. 
defines the embedded hypersurface E e ). In terms of coordinates, this gauge transformation 
corresponds to e-dependent coordinate changes y l (y\e) of the intrinsic coordinates in the 
hypersurface. It is obvious that this gauge transformation does not affect the first and second 
order perturbed metrics K', K". However it does affect how points with fixed coordinates 
y l move on spacetime and therefore it affects the vectors Z\ and Z 2 - In order to find their 
gauge transformations, we need to define suitable first and second gauge vectors u~[ and u 2 . 
It is clear from our previous discussion that they can be defined as 



Ui 



dXe 



de 



«2 



9 2 (Xe+h O Xe *) 



e=0 



dhde 



=7i=0 



Being defined on E they can be pushed-forward to spacetimes vectors defined on E . As 
before, we shall use the same symbol for both objects. 

Proposition 6 Under a gauge transformation on E defined by gauge vectors u\ and u 2 , 
Z\ (p) and Z2 (p) at any point p G Eq transform as 



#i(p) +«i(p), 



z 2 (9, ( P ) 



Z 2 + u 2 + 2VuiZx - afaijuiur 1 



11 



(48) 
(49) 



where n is a unit normal, (n, n) = a, and Kij is the second fundamental form 0/E0. 

Remark: The covariant derivative in the second expression is a spacetime covariant deriva- 
tive, not a covariant derivative on E push-forwarded to £0 

Proof: A similar calculation as the one leading to (J47j) shows that 



d 2 x\{y) 



de 2 



+ rf\yW(yW(y), 



(50) 



where T^* are the Christoffel symbols of (E, h). Writing (jllj) with the substitutions $ e 



[9) 



*e O Xe, Zi 



Zi {9)a (y) 



immediately implies 



de 



dx l e(y) 



e=0 



Oe 
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where ef(y) = d ^ d ^ are tangent vectors to S and we have used that Xe=o = Is- This 
proves ([48)1 . For the second derivative we find, after gauge transforming (|12p. 



e=0 



+ 

2/ 



+ + r ^^(^/))efej) • (51) 

The term in parenthesis can be rewritten as efV^e". Decomposing this vector into tangent 

and normal components we have efV^e" = —aKijn a + ef, where we have used the 
fact that the projection on So of the spacetime covariant derivative is precisely the three- 
dimensional covariant derivative. Using this and (|50j) into (|51[) yields the result. □ 
Remark: The gauge transformation on E only affects the way how points on £ are 
identified before the e-derivative is performed. One could think that the only effect of this 
should be transforming the vectors Z\ and Z 2 with tangential components. While this is 

clearly so for the first order vector, it is not true for the second variation (i.e. Z 2 9 — 
Z2 contains in general normal components). The reason is that the vector Z2 measures 
essentially spacetime accelerations, and curves fully contained within So in general have a 
non-zero normal acceleration unless the hypersurface is totally geodesic, i.e. = 0. Notice 
that in this case the transformation (J49|) does not have normal components (provided Z\ is 
tangent to S , of course). 



8 First and second order linearized matching condi- 
tions 

In this section we apply the previous results to perturbed matching theory. More specifically, 
we will consider two spacetimes (Ai + ,g + ) and , g~) joined across a common hypersur- 
face S , i.e. such that the so-called matching conditions are satisfied on S . On both regions 
we perturb the metric to first and second order with K ;± and K n± respectively. Our aim 
is to obtain necessary and sufficient conditions on So such that the matching conditions are 
also satisfied in a perturbed sense (i.e. to first and second order). 

Let us start with a brief discussion of the matching conditions. Matching theory between 
spacetimes deals with two C 2 spacetimes (A4 + ,g ± ) with boundary. The respective bound- 
aries are C 3 hypersurfaces of M.^- which are called matching hyper surf aces. We will denote 
them by £ . Although a fully successful theory can be developed for boundary of an arbi- 
trary causal character [3U| , for simplicity we will concentrate here in the case where both E ± 
are either timelike or spacelike everywhere. The matching theory asserts that a C° spacetime 
(M. 1 , g l ) can be constructed by joining (A4 + , g + ) and (A4~ , g~) if and only if there exists 
a diffeomorphism (p : S + — ► £~ which is an isometry with respect to the induced metrics 
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/i ± on each boundary. This is equivalent to introducing an abstract (n — l)-dimensional C 3 
manifold E, and demanding the existence of two embeddings <I >± : E — > M.^ such that (i) 
$ ± (E) = C M.^ and (ii) that the two induced metrics /i ± = $ ± *(5' ± ) on E coincide, i.e. 
h + = hr . Condition (ii) is called first set of matching conditions. Furthermore, the Riemann 
tensor in the joined spacetime (.M*, #*) is free of Dirac delta distributions if and only if the 
second fundamental forms on E coincide, i.e. 

$+* (V+n+) = $-* (V _ n _ ) , (52) 

where ± always means that the objects are calculated from the A4 ± side. The unit normal 
to E 1 * 1 are chosen to have the same relative orientation after the matching, i.e. either n + 
points outwards and n~ inwards or vice versa ("inwards" and "outwards" are well-defined 
concepts for vector fields on the boundary of a manifold-with-boundary). Notice that the 
matching conditions are fully covariant both with respect to the spacetimes Ai ± and with 
respect to the matching hypersurface E. This means that any local coordinate system in 
M. + and any coordinate system in Ai~ are equally valid to impose the matching conditions. 
Moreover any coordinate system can be used in the abstract manifold E. 

We want to perturb the spacetimes with boundary (M. ± 1 g ± ) and analyze under which 
conditions the matching conditions are satisfied perturbatively (assuming the background 
spacetimes do match across their boundary). First of all, we need to consider how to define 
perturbations for manifolds with boundary. As mentioned above, in important ingredient in 
perturbation theory is the need to identify the different manifolds M. t with each other. Now, 
our manifolds M.+ have boundary (we concentrate on the "+" side; similar considerations 
hold for M.~). If we imagine them as subsets of larger manifolds without boundary M e and 
identify these, it is clear that, generically, the identification will not transform the boundaries 
E e among themselves. We could insist that the identification preserves the boundaries, but 
only at the expense of restricting strongly the gauge freedom (at least near the boundary). 
This may be quite inconvenient for other purposes. Thus, it is better to let the boundary 
"move" freely in the identification. From the point of view of the family of manifolds with 
boundary this means that, strictly speaking, we are not taking diffeomorphisms between 
them. They are diffeomorphisms except in some closed neighbourhood of the boundary 
E+, for each e. Thus, in strict terms, we cannot talk of a background manifold M.$ with 
boundary E^ and a family of metrics g e defined on it. Nevertheless, in perturbation theory 
we only care about derivatives at e = of the e-family and this can be consistently defined 
on points at the boundary Eq by taking one-sided derivatives (i.e. restricting the variations 
to positive or to negative e depending on the point of the boundary we are considering). 
This allows one to define, similarly as in the case without boundary, a background manifold 
with boundary M. + and two symmetric tensors K' + and K" + defined everywhere up to 
and including the boundary, so that we have a proper perturbation theory up to second 
order (or higher order if desired). Having this in mind, we will abuse notation and still talk 
about diffeomorphisms between different M.^ . This also allows us to talk about how the 
hypersurfaces E e move on the background and therefore introduce vectors Zi + and Z 2 + on 
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the unperturbed boundary, exactly as we did in Section HJ 

With this particularity in mind, it now easy to write down the perturbed matching 
conditions. Indeed, for each e the matching conditions demand the equality of the first 
and second fundamental forms on each side, i.e. hf = h~ , = k~. This tensors are all 
defined on the abstract hypersurface S and therefore can be compared with each other (and 
differentiated with respect to e). It is also clear that the matching conditions will be satisfied 
in a perturbed sense if and only if the first and second derivatives of the first and second 
fundamental forms coincide from both sides. Using the explicit form of these derivatives 
found in previous sections we have a practical method of determining whether perturbations 
of a spacetime constructed by joining two regions across S can be matched across this 
hypersurface. We state this result in the form of a theorem 

Theorem 1 Let (A4,g) be a spacetime constructed by joining two spacetimes with boundary 
(Ai + ,g + ) and (A4~,g~) across their corresponding boundaries £ . Let £ be an abstract 
copy of S + and $ ± : S — > A^ ± be the embeddings defining the background matching. Let 
also K'^ and K" be first and second order metric perturbations in M. ± . 

The first order perturbed (i.e. linearized) matching conditions are fulfilled if and only if 
there exist two scalars Qf and two vectors T\ on £ for which 

h' + — h'~ k' + — k'~ 

ij iji ij ij ' 

holds, where h!^ , k'^ are given in Proposition^ after the substitution Q\ — > Qf, T\—>T\ , 
g -> g ± , K' -> K r± and ef -> e?* 

The second order perturbed matching conditions are satisfied if and only if there exist two 



■ ± 



scalars Qf and two vector fields T 2 on £ such that 

ij ij' 1 ij ij' 1 

where these objects are obtained from ^U^ -^l]) after similar substitutions. 

Remark: It is important to stress the fact that satisfying the perturbed matching conditions 
require the existence of vector fields Z\ and Z 2 such the equations above are satisfied. 
These vector fields are not known a priori. Moreover they need not be the same vector on 
both sides. This is obvious from the fact that these vectors are gauge dependent and the 
gauge may be chosen differently in the different regions Ai (actually one can not even 
compare the two gauges, in general). The gauge can always be chosen so that these vectors 
coincide but this may not be the most convenient choice. Linearized matching conditions 
have often been analyzed by using specific gauges where the vectors Z\ and Z 2 take simpler 
forms. A common choice is to use Gauss coordinates adapted to the matching hypersurfaces 
for all e (this obviously makes Z\ = Z\ = 0) and then transform to the desired gauge. This 
is the approach taken in jSH] for instance. In spherical symmetry, the linearized matching 
conditions in arbitrary gauge was first studied in |37j, [HH] and completed in The 
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general linearized matching conditions in an arbitrary background were first presented by 
Mukohyama [22] • in this paper a vector field Z was introduced describing the perturbation 
of the matching hypersurface to first order. However, this vector was assumed to be the 
same in both sides M. + and M.~ . As I have already stressed this need not be case and 
the fully general perturbed matching conditions require the use of two vectors Z\ and Z\ 
(and two more vectors to second order). Notice also that the gauge freedom within £ (see 
Proposition IHJ) allows us to choose the tangential components of Z\ and Z 2 in any way we 
want, but only on one of the sides, i.e. either on Ai + or on . Once a choice on one side 
has been made, the other side must be left free and determined by the matching conditions 
(if they happen to be consistent). This is similar to the fact that when solving a matching 
problem one not only looks for a pair of matching hypersurfaces with suitable properties, 
but also for a specific pair of embeddings on each side, i.e. a way of identifying the two 
hypersurfaces pointwise. 

In this theorem, only non-null hypersurfaces are considered. This is because the classical 
Darmois matching conditions (discussed above) are not adequate for hypersurfaces with 
null points. In that case the continuity of the second fundamental form does not ensure 
the absence of distributional parts in the Riemann tensor. The matching conditions for 
null hypersurfaces where first discussed by Clarke and Dray jHE] and later extended to 
hypersurfaces of arbitrary causal character (including a changing one) in jHO] • They involve 
the continuity of a tensor that generalizes the second fundamental form. In order to find 
the perturbed matching conditions in this case we would need to find how this new tensor 
is perturbed to second order. This issue is of interest and should be studied. The methods 
described in the present paper are useful to find perturbations of any geometric tensor defined 
on a hypersurface and therefore are applicable to this situation too. The calculations for 
hypersurfaces with null points are probably more difficult but still manageable. 
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9 Appendix A: Proof of Proposition El 

Proof: Let us start with h" . From Lemma Q we find 

d e d e h t \ e=0 = $* [K" + 2£ Z ~K' + £^g + £^£^g) . (53) 
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Applying (J53j) with B — > g and using (p3|) with Fi = F 2 = Q\ we readily obtain 

-2T 1 (g 1 )) K a(i + 2Q? (-nWRatf,, + K^ng) + 2aD a Q 1 DpQ 1 + n a P^ + n^P a . 

Here and in the following P a , P a p, • • ■ stands for expressions whose explicit form does not 
concern us. Notice that its meaning may be different even in different parts of the same 
formula. 

Substituting into we observe that a term 2£f (K' + £ g g) appears. From Proposition 
121 the pull-back of this term is just 2£fh'. Recalling that Z 2 = W + V g Z\ and the 
decomposition Z 2 = Q 2 n + T 2 yields the first two terms in (jjOJ). Only 2£q x ^K' ' a p remains 
to be analyzed. This is dealt with using Lemma El with X — > Q\n which gives 

2£ Qlfi K' a p = -AQxn^S'^p + 2£ Q ^,g a p + AoQ x Y' ' K al3 + n a P p + npP a , 

and expression (|4Hj) follows directly. Let us next consider k", which involves the longest and 
most difficult calculation. Applying Lemma Q to n e we find 

2d e d e K e \ e=0 = (£ n2 g + 2£ Hl K' + £ n K" + 2£ fi £ Hl g + 2£ fi £ H K'+ 

+ £ £ fig + £ g x £ z% £n9) ■ (54) 

In Proposition |21 we evaluated d e K e \ e= o, which required calculating the pull-back on X of 
M = £f i £^g + £ft 1 g + £nK' (see (j28j) ). We want to identify in (jHljl terms giving the Lie 
derivative of M along T\. Adding and subtracting 2£f £g £^g we can write 

2d e d e K e \ t=Q = <T (£ H2 g + 2£ Hl K' + £ n K" + 2£ QlH £ Hl g + 2£ QlH £ H K'+ 

~\~ £ \y £ ag £ Z\ ^ Z\ ^ ™9 2£rp^£ g^£fig -\- 2<£j?ikf) . (55) 

Let us deal with the different terms in this expression, starting with £ J p L K'. A direct appli- 
cation of P0|) and the explicit expression for n\ (j^H) yield, after using 

£n 1 K' a p = -Y'\ a/3 - °-Y'£ ? g afi + - °- (r' a D^Y' + r> ' pD a Y') - 

- £ ?+Q 1 3+a^Q 1 K '^ + n «Pp + TlpP a . (56) 

Next, we analize the combination 2£q 1 ^£ 7 ^ l g + 2£q x ^£^K' . Using expression (J31|) and the 
first equality in we get 

2£Q 1 ft£ri 1 9a/3 + 2£q^£^K' a p = (7 £ Q 1 rf£ Y'n9a/3 ~ 2 £ Q 1 fi£ Q 1 a+a S r a dQ9af3 ~ ^£q 1 u { n fi^'ap) • 
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Using now identity (jzoj) with F± = Qi, F 2 = K' and writing the Lie derivative of the last 
term explicitly using covariant derivatives yields 

^£Q 1 n£n 1 9al3 + 2£q^£^K' 'a/3 = < J £Q 1 Y'a9af3 

+2aQ 1 Y' (-nWRcMv + «a^«/) + & a Qf7 p Y' + V a Y'V >Qi) + (57) 
+2o-K Q/3 Qira(F / ) - 2£ Qlil £Q lS+r7gr! , d Q 1 g a { 3 - 4Q 1 a M <S / ^ - 4Q 1 n /i n I/ V 1/ «S / ^ - 
-471X5^ V/»Qi - 4nX^V a g! - 4n M <S'^^ - 4^5'^^ + n Q P^ + n p P a . 

Next we analyze the terms involving third derivatives in (|55|). Lemma |S] implies 

£ W £ n9 £ Z\ £ Z\ £ft9 '^'£'Ti £ Z\ £ft9 
Z2 ™9 £'f 1 £j ! 1 £n9 C1H+D1 £ft9 £Q\a£n9 £ Q\fi£ Q 1 n£fi9 ■ 

Except for the last term £Q 1 n£Q 1 n£fig, the pull-back of the right-hand side is easily obtained 
from Lemma El and ()23|) . For the last term we have, directly from 



£Q 1 n£Qrn£ri9 — £qiH£qiS9 + 2Qifl(Qi) (~ n^n" R afl f3 u + K^K^) + 2Q\£q^ {K all Kp — 

n^n v R a ^ u ) + 2(7(5! ( a M «a M {Qiap + oDpQ x ) + a^K^ (Qia Q + crD a (5i)J + ™ a P/3 + ^^.(58) 

In order to elaborate this expression further we expand the Lie derivative in the third sum- 
mand in terms of covariant derivatives. We get 

n^rf £ QlH R a ^ u = Qin^n u n s V &R a ^ v + 2n(Q l )n fJj n l, R afl[3l/ + 
Qi (a fl n u + Ka-Rff^n^n" + K^Rg^nV) + n a P p + n p P a . (59) 

For the Lie derivative of n a p along Q\n we have 

£Q 1 nK-a(3 = -crQia a ap + Qi {-n^n v R ail p v + K afl K^) + Q-h^h^zg^, (60) 

which follows from £Q lf thp = —aiQiap + aDpQ^rf 1 and K Q/ g = ^h^h%£ 'q^Q '\w after applying 
((211). Inserting (JHH)) and (jHOj) into (JH%|) the following expression is found 

£Q 1 n£Q 1 n£n9af3 = £ Q x n£ Q^a9af3 + 2Ql^(Ql) ( — U^Vb' R ail ^ v + K^K^) 

-2Q2 (^n"n 5 VA^, + a»n v R ail p v + nWR^ + 2^R s ^n^n u + 2k$ R s ^n v ) 

+2Q 1 a fl K flP D a Q 1 + 2Q 1 a fl K fia DpQ 1 + Q\ (n^£ d g afl + K££ s gp fl ) + n a P p + n p P a . (61) 

It only remains to calculate the first term in (|55|) . i.e. to find the vector n 2 = d e d e n t \ t= Q. 
The calculation is somewhat long and will be given in Appendix B. The result is 

n 2 a = n a (~<rY" + ^Y' 2 + ar' - r" a - (Q 2 a a + aD a Q 2 ) + aFV" 

+ (da a + aD a d) - (n(Q{) + oY') (Q ia a + aD a Q 1 ) + 2K' taP (r'p + Q^p + aD p Q 1 ) 

+ [Qi^Qia + agradg^ + 2Q?a^ + 2aQ 1 K a tl D"Q 1 . (62) 
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Our interest is to calculate £n 2 9p,u- Only the term involving 2K' ta ^ \t' ' p + Q±ap + aDpQi) 
requires further analysis. Applying Lemma El and the fact that S'^p = ^{K' 1 )^ + o-T^K a p + 
'Pap+naPp* 1 +npP a fJ- , which follows directly from the definition of S' and the decomposition 
flUD , we get 

£2K' t ^{T'„+Q 1 a„+aD v Q 1 )9al3 = ^r'+QiS+crgi-adQi K> + 

+4 (t'^ + Qia^ + aD^Qx) (S'^p - ^K a p) + n a Pp + npP a . 

We are now in a position where all terms in (|HH)l can be collected. Several obvious can- 
cellations happen which will not be described in any detail. More subtle is the use of the 
following identity 

~£'Qxa,£'Qxn9aP ~ £Q\a,£n9aP + Q\^£a9pp + Ql K p£a9ap + %Ql K^pD a Qi + d^K^DpQx) 

-2Ql {n^a v R ail p v + a^n u R a ^p u ) + 2Q 1 a(Qi)K a p + 2£ Q 2 a u K ^g a p = n a Pp + npP a . 

This expression follows by direct calculation using the Codazzi identity written in the space- 
time form D a K^ u — D^K av = n s R^spph^h^h^ and the fact that D a ap — Dpa a = 0, which is a 
direct consequence of the definition of acceleration in our hypersurface orthogonal case. 

Finally, in order to arrive at the final expression, two more ingredients are required. The 
first one is 

[Q^gradQi] = Q 1 D a {n{Q 1 )) - <rQ 1 n{Q 1 )a a - 2Q 1 k^D"Q 1 
- (D^QxD^Qx + crQia(Qi)) n a , 

which is checked directly. The second one is n(Y') = 2r' a a a + 2n fJi n p n s S'p S , which is imme- 
diate. Using also the explicit expression for C\ and D\ in (J35J) and collecting all terms we 
find dHJ. □ 

10 Appendix B: Calculation of n<i 

In this appendix we find an explicit expression for d e d e n t . 

Lemma 6 With the same notation and conventions as in Proposition we have 

n 2 a = n a i-\oY" + \y' 2 + ar'rA - r" a - (Q 2 a a + aD a Q 2 ) + aY'r' a 



+ (C ia a + aD a C l ) - (n(g0 + o-Y') (Q ia a + oD a Qi) + 2K' taP (r'p + Q x ap + aDpQ x ) 



Q 1 n,Q 1 a + agr&dQ 1 + 2Q?a"< + 2aQ 1 K™D tl Q l . (63) 
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Proof: We proceed as we did for rfi, i.e. we first determine the normal component of ri2 and 
then its tangential part. Taking the second e derivative of (n t ,n t ) ge = a and evaluating at 
e = we find 

2n 2 % + 2nSm„ + An^rfK'^ + n»n v K\ v = 0, 
which after substitution of the expression for n\ (}3T| yields 



1 



n 2 "n IM = --Y" + -oY' 1 + (r^ - - aD,Q{) (r' M + + aD»Q x ) . 



(64) 



For the tangential components, it is convenient to use the second variation of the normal 
one-form, i.e m.2 = d e d e n e \ e= Q. The relationship with is immediate from the second 
e-derivative of g e (n e , •) = n t , i.e. 

n 2 a = -K" aP n p + aY'K ,aP n p + 2K ,a<i (r 1 'p + Q x a p + OpQx) + m 2 a . 

Decomposing this into tangential and normal components and using one finds 



n 2 a = n a 



(65) 



(66) 



-r + oY'r' a + 2K'^ (r'p + Q x a p + aD^) + h%w^ . 
It only remains to find h^m^ ■ Applying Lemma ^ to Q e (n e ) = 0, Ve yields 

{£yv n + ■^z 1 ^z 1 n + r 2£'z 1 ra ^ + ra ^) = 0- 
Moreover, Lemma applied to n leads to 

~£Q\a n + £Qin£Q r n n - 

Now, for any pair of functions F\ and F 2 we have £ {F2T1) = n(FiF2)n a + F2{F\a a + 
aD a Fi). Using also the fact that §*(£^n) = for any vector field V tangent to E, i.e. 
£^n oc n one finds 

£w~n a + £^£^n a = Q 2 a a + aD a Q 2 + 2£f x (Q\a a + aD a Qi) - (Cia a + aD a Ci) + 

{Q\a a + aD a Qx) + £ Ql , H (Qia a + (?D a Qi) + Pn a . 

Moreover, from (JHUj) 



£ z ~ 1 m ia = ~£f 1 (Qi a a + vD a Qi) + -1" + o-DaQx) - £ QlH (Qia a + aD a Qi) + Pn a . 
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Then (|5l)j) implies 

him 2 p = - (Q 2 a a + aD a Q 2 ) + (da a + aD a d) - (n(Q x ) + aY') {Q x a a + o-D a Qi) + 

+hP£ QlH (Q ia p + aDpQr) . (67) 

Let us elaborate the last term. Using £ Qlf ih al3 = —2QiK af5 — an a (Qia f} + aD^Qi) — 
an^{Q\a a + aD a Qi) 1 and integrating by parts we find 

h a ?£ QlH (Qtap + aDpQi) = [Q x n, Q x a + agrad Q 1 ] a + 2Q 2 1 a^ a fM + 2aQ 1 n a ^Q 1 

+an a (Q lflA1 + aDM (Qia" + aD^Qt) . 

Plugging this into (jHTjl and the resulting expression in (|fi5j) . the Lemma follows. □. 
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